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ABSTRACT 


i  I 


An  approximate  solution  is  presented  for  determining  transient 
temperature  distributions  and  ablation  rates  for  solid  spheres  of  low  thermal 
conductivity  and  constant  thermal  properties,  subjected  to  point  symmetric 
aerodynamic  heat  inputs.  A  short  time  solution  for  the  temperature  which 
is  valid  prior  to  ablation,  is  developed  by  use  of  Laplace  trans'i.  rms .  The 
ablation  solution  is  derived  by  approximation  of  the  radial  temperature  pro¬ 
file  by  an  exponential  function  and  utilization  of  the  heat  balance  technique. 
This  results  in  a  second-order,  nonlinear,  L-'dinary  differential  equation 
for  the  position  of  the  ablating  surface. 

Results  of  the  numerical  integration  of  the  approximate  solutions 
are  compared  with  experimental  data  obtained  at  the  stagnation  point  of  a 
Plexiglas  nemisphere  tested  under  hypersonic  flow  conditions.  A  discussion 
of  the  techniques  employed,  as  well  as  the  comparison  with  experimental  data, 
appear  in  the  body  of  this  report. 
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SYMBOLS 


initial  radius  of  sphere 

speci f ic  heat 

heat  transfer  coefficient 

thermal  conductivity 

heat  of  sublimation  or  melting 

stagnation  pressure 

constant  heat  flux 

radial  coordinate 

radius  of  ablating  heated  surface 

nondimensional  radial  coordinate,  r/a 

time 

melt  time 

2 

nondimensional  time,  nt/a 
temperature 

heated  surface  temperature 

initial  temperature 

melt  or  sublimation  temperature 

radial  temperature  profile  from  pre-melt  analysis 

stagnation  temperature  of  flow  field 

nondimensional  temperature,  (T  -  T.)/(T  -  T.) 

transformed  variable,  rT  or  r(Ts  -  T) 

Laplace  variable 


*1 


Ii 


nond  i  mens  i  ona  I  heat  transfer  parameter,  aQ.Q/ pin 


thermal  diffusivity 

nond i mens i ona I  ablation  radius,  r  /a 

s 

nond i mens i oim  1  material  physical  parameter, 
dens i ty  of  solid 


n(t  -  t  ) 


nond imens i onal  time,  - 2 - 

a 

nond imens i ona 1  heat  transfer  parameter, 


ah 


MT  '  T; ) 

m _ i _ 

pin 


indicates  differentiation  with  respect 


to  T 


INTRODUCTION 


The  high  inputs  and  aerodynamic  forces  encountered  by  a  space 
vehicle  upon  re-entry  into  the  atmosphere  have  created  numerous  areas  of 
interest  for  the  scientist  and  engineer.  One  area  of  recent  emphasis  has 
been  ablation,  or,  the  utilization  of  the  latent  heat  of  vaporization  or 
sublimation  of  the  material  in  order  to  reduce  the  heat  input  to  the 
vehicle  structure. 

Considerable  research  has  been  performed  in  problems  of  transient 
temperature  distributions  in  bodies  undergoing  phase  changes.  References 
[1]  through  [7],  together  with  the  lists  of  references  therein,  contain  a 
fairly  comprehensive  bibliography  of  this  subject. 

The  current  research  was  motivated  by  the  desire  to  obtain  a  pre¬ 
diction  of  ablation  rates  and  temperature  profiles  for  thick-walled  hemispheres 
of  low  thermal  conductivity  subjected  to  aerodynamic  heat  inputs.  Tests  of 
Plexiglas  and  Nylon  hemispherical  models  are  now  being  conducted  in  the  hyper¬ 
sonic  tunnel  of  the  Polytechnic  Institute  of  Brooklyn.  The  experimental 
results  available  on  one  of  the  models  are  compared  to  the  present  theory. 

A  more  extensive  comparison  will  be  the  subject  of  a  future  report. 

5n  this  report  an  approximate  theoretical  method  is  developed  for 
determining  transient  temperature  di stributions,  together  with  ablation  rates 
for  solid  spheres  of  low  thermal  conductivity  and  constant  thermal  properties 
subjected  to  point  symmetric  aerodynamic  heat  Inputs.  Since  the  high  thermal 
Input  conditions  encountered  In  re-entry  produce  steep  radial  temperature  gra¬ 
dients  at  the  heated  surface  of  a  material  of  low  thermal  conductivity,  the 


2. 


material  temperature,  at  a  short  distance  from  the  heated  surface  is  essen¬ 
tially  the  same  as  the  initial  material  tern: e rat u re.  Consequently,  the  solu¬ 
tion  for  a  spherical  shell  of  sufficient  thickness  may  be  approximated  by 
the  solution  for  a  solid  sphere.  The  problem  as  developed  in  this  report 
is  divided  into  two  parts:  (1)  the  pre-melt  solution,  and  (2)  the  ablation 
solution. 

The  pre-melt  solution  utilizes  Laplace  transformations  which  re¬ 
sult  in  a  short-time  solution  to  the  heat  conduction  equation  [8].  This 
solution  is  considered  valid  until  the  sublimation  or  melting  temperature 
of  the  body  is  reached  on  the  heated  surface.  The  determination  of  the 
time  at  which  the  above  condition  occurs  is  of  primary  importance  in 
correlating  experimental  ablation  data  with  the  ablation  theory. 

The  ablation  solution  is  approximated  by  employing  the  Goodman 
heat  balance  technique  [2].  An  exponential  temperature  profile  is  assumed 
in  the  unablated  sphere  and  substituted  into  the  heat  conduction  equation. 

The  resulting  equation  is  integrated  over  the  region  of  the  solid,  and 
yields  ar  ordinary,  nonlinear,  second-order,  differential  equation  relating 
the  ablation  depth  to  time.  The  latter  equation  is  solved  by  numerical 
techniques. 

The  ablation  solution,  used  together  with  the  pre-melt  solution, 
provides  a  simple  and  reasonable  approximation  for  the  ablation  of  low 
conductivity  spherical  shells. 
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PRE-MELT  SOLUTION 


A  short-time  solution  for  the  temperature  distribution  in  a  solid 
sphere  subjected  to  a  point  symmetric  aerodynamic  heat  input  is  obtained 
by  utilizing  Laplace  transformations.  The  appropriate  form  of  the  heat 
conduction  equation  is 


dT 

at 


2  aT  v 
7  ar 


with  the  conditions 


0  <  r  <  a 


(0 


h(Ts  -  Ta)  -  k(3T/3r)r  .  , 

T(r,0)  -  T. 

T(0,t)  is  finite 

where  the  sumbols  are  defined  as 

a  initial  radius  of  sphere 

h  heat  transfer  coefficient 

k  thermal  conductivity 

r  radial  coordinate 


(2a) 

(2b) 

(2c) 


t 


time 


4. 


T 

T 

a 

T. 

i 

T 

s 

K 


temperature 

heated  surface  temperature 
initial  temperature  of  material 
stagnation  temperature  of  flow  field 
thermal  diffusivity 


Equations 
the  substitution 


(I)  can  be  transformed  into  a  more  convenient  form  by 


u  =  r(T$  -  T) 


The  result  is 


and  equations  (2)  become 


0  <  r  <  a 


(3) 


du 

dr 


=  0 


at  r  =  a 


(4a) 


u(r,0)  =  r(Ts  -  Tj )  (4b) 

u(0,t)  =  0  (4c) 


where  ^  =  (ah  -  k)/k. 
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The  Laplace  transformation  is  now  applied  to  Eq.  (3) 


J  .->*  U  at 

o  3r 


and  by  using  Eq.  (4b),  the  resulting  subsidiary  equation  is 


<Ts  -  Tl>r 


where 


q  =  p/H  . 


Equations  (4a)  and  (4c)  become 


dr  a 


at  r  = 


u  =  0 


at  r  = 


A  general  solution  to  Eq.  (5)  is 


(T  -  T.)r 

u  .  C,eqr  .  C,.-^  +  - 


where  Cj  and  Cj  are  constants.  These  constants  are  evaluated  by  applying 


6. 


Eqs.  (6a)  and  (6b).  The  result  in  series  form  is  as  follows: 


u  = 


(VT,)r 


0+l»)(Ts-T|) 


y  (  ))"  (q-u/a)n  f  -q[(2»»l)a-rl  e-q[ (2n+l)a+r)\ 
L  (qVa)n+l  1  ; 


n*0 


(7) 


For  a  short-time  solution,  Eq.  (7)  may  be  truncated  after  the  n  =  0  term. 

For  the  range  of  parameters  considered  in  this  report  the  higher  order  terms 
oc  Eq.  (7)  can  be  shown  to  be  negligible. 

The  inverse  transformation  of  Eq.  (7)  for  n  =  0  is 


t' 


Kt 

2 


r'  =  r/a 


where 
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Equation  (8)  represents  the  short-time  solution  for  the  temperature-time 
history  at  any  point  on  or  within  the  sphere. 

In  order  to  correlate  the  experimental  ablation  data  with  the 
ablation  theory,  the  time  (tm)  at  which  the  melting  temperature  (Tm)  of 
the  material  is  initially  reached  on  the  heated  surface  must  be  deter¬ 
mined.  Surface  temperatures  versus  time  for  various  values  of  the  heat  input 
parameter  (p)  have  been  calculated,  using  Eq.  (8),  and  the  results  are 
shown  in  Fig.  I.  For  a  given  material  and  heat  input  condition  Fig.  1 
can  be  used  to  determine  the  melting  time  (tm)> 

Once  t  has  been  determined,  the  radial  temperature  profile  at 
t  *  t  can  be  calculated  from  Eq.  (8).  The  determination  of  this  profile 
is  not  essential  for  the  solution  of  the  ablation  problem,  but  it  must  be 
shown  here  as  a  justification  for  the  simplified  profile  chosen  for  the 
ablation  theory.  One  such  profile  is  presented  in  Fig.  2  as  a  representative 
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ABLATION  SOLUTION 


The  transient  temperature  distributions  and  ablation  rates  for 
the  solid  sphere  are  approximated  by  employing  the  heat  balance  technique 
due  to  Goodman  [2].  The  solution  is  considered  valid  from  the  time 
melting  commences  on  the  heated  surface  to  the  time  that  the  predicted 
temperature  increase  at  the  center  of  the  sphere  ceases  to  be  a  negligible 
percentage  of  (T  -  T.).  The  initial  assumptions  made  in  this  analysis 
are  as  follows: 

a.  the  heated  surface  remains  at  the  melt  or  sublimation 
temperature  (T^). 

b.  the  melt  or  products  of  sublimation  are  immediately  removed 
upon  formation, 

c.  the  thermal  properties  of  the  solid  are  independent  of 
temperature. 

d.  the  heat  flux  (Q.o)rematns  constant1. 

The  appropriate  form  of  the  heat  conduction  equation  is 


*1  _  tA  .  2  <yn  .  0  <  r  <  r  (t) 
dt'V  r3rJ>  t>t 


(9) 


The  function  rs (t )  denotes  the  time-dependent  radial  coordinate 
of  the  receding  outer  surface  of  the  sphere. 


In  the  case  of  aerodynamic  heating  this  assumption,  together  with  the  assump¬ 
tion  that  the  heated  surface  remains  at  a  constant  temperature,  is  analogous 
to  assuming  that  the  heat  transfer  coefficient  (h)  is  constant. 
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The  boundary  and  initial  conditions  to  be  satisfied  are 


T(r  ,t)  =  T 
s  m 

(I0») 

T(r,t  )  =  T  (r,  t  ) 
m  o  m 

(10b) 

T(0, t)  is  finite 

(!0c) 

r5(tn,>  “  3 

(I0d) 

-"l3T 

r  ,t 

(10c) 

where  TQ(r,tm)  the  temperature  profile  at  time  (t^)  determined  from  the 
pre-melt  analysis,  and  p  and  L  are  the  density  and  heat  of  sublimation  of 
the  material,  respectively. 

Equation  (9)  is  converted  to  a  more  convenient  form  by  applying 
the  transformation 

u  =  rT 

whereupon 

.2 

■*1f>  °<  >■<  v  t>  t  .  (ii) 

9r 

I 


^u 

at 


tit 
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integrating  Eq.  (11)  over  the  volume  of  the  solid, 


or 


3rJ  S  udr  -  “(vl)  3T‘  K[(lr) 


V‘ 


0  1 

0,t 


(12) 


Reintroducing  Tr  for  u,  Eq.  (12)  may  be  expressed  in  the  form 


£  J0S  Trdr  ■  rsT(vt)  d r =  x  f  +  T(rs-t)  -  T(0>t)] 


(13) 


Condition  (10b)  is  now  satisfied  approximately  by  assuming  the 
temperature  profile 


T  -  T. 

i 

T  -  T. 

m  i 


r  rQ.  +  pL(dr  /dt)-,  . 

■  “*•  {{  °  k(-T  ztf— ](r,  -  D} 


(14) 


m  i 


This  temperature  profile  also  satisfies  conditions  (10a),  (10c)  and  (lOe). 
Substitution  of  (14)  into  (13)  yields  the  following  equation: 


[2(1  -  e~or)  -  a(l  *  e-®)] 


d2r 


dt 


r  dr 


■  £  <r>3<Tm  -  Ti>K« ♦  c  -  e'°)] *  [•  -  C  *  •->]  (f>d f}  <l5> 


11. 


where 


12. 


_  f^> 

^  ”  Plx 

k<Tm  -  V 

^  pLx 


H  = 


£D 

dT 


n  = 


A  stepwise  numerical  integration  procedure  is  used  for  the  solution  of 

Eq.  (17);  starting  with  the  initial  conditions  on  1)  and  T)  .  It  can  be 

shown  from  condition  (lOd)  that  7]  at  t=0  is  unity,  and  for  continuity  of 

the  heat  flux  at  the  so-called  melt  time,  that  7)  at  t  =  0  is  zero.  These 

•  •  • 

initial  values  of  7]  and  7]  are  now  used  to  find  7)  (0)  from  Eq.  (17) •  A 
value  of  the  time  interval  (At  =  t2  -  t^)  is  chosen,  and  the  following 
Taylor  series  expansions  are  used  to  proceed: 


2 

T)(t2)  =  TJ(t j )  +  T)(t , ) At  h-’tKt,) 

T)(t2)  =  TKt,)  +  1)’  (t,)At 


(18) 


The  values  of  7](t2)  and  T|(T12)  thereby  obtained  are  substituted 
•  • 

into  Eq.  (17)  to  obtain  7](t2).  This  process  is  then  repeated  for  each  suc¬ 
cessive  time  interval. 

The  resulting  time  history  of  the  ablating  surface  (T)  vs  t)  for 
various  values  of  the  heat  input  parameter  (y),  and  two  values  of  the 
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material  parameter  (cp),  are  shown  in  Figs.  3  and  4.  The  corresponding  time 
histories  of  the  ablation  rates  Cty  vs  t)  are  shown  in  Figs.  5  and  6.  The 
values  of  <p  used  correspond  to  the  physical  properties  of  Plexiglas 
(<p  *  0.511)  and  Nylon  (<p  =  0.108). 

The  values  of  the  ablation  rate  and  ablation  radius  so  deter¬ 
mined  can  be  substituted  into  Eq.  (14)  to  yield  the  temperature-time 
history  at  any  point  within  the  body. 
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COMPARISON  WITH  EXPERIMENT 

The  theoretical  results  of  the  preceding  sections  are  now  com¬ 
pared  with  an  experimental  result  obtained  at  the  stagnation  point  of  a  Plexi¬ 
glas  hemisphere  tested  in  the  hypersonic  tunnel  of  the  Polytechnic  Institute 
of  Brooklyn.  The  shroud  technique  was  utilized  to  produce  pressure  distri¬ 
butions  and  heat  transfer  rates  corresponding  to  high  speed  re-entry  of  the 
hemi sphere. 


Dimensions  of  Model 


Outer  Radius,  a 

3.875  in 

Inner  Radius,  b 

2.875  in 

Test  Conditions 

Stagnation  Pressure,  Ps 

19.5  psia 

Stagnation  Temperature,  Ts 

0 

0 

0 

30 

Heat  Transfer  Coefficient  at  Stagnation  Point,  h 

0.0,55  2BTUo 

Ft  sec  R 

Initial  Temperature  of  Model,  T. 

540°R 

Duration  of  Test  Run,  t^ 

37  sec 

Ablation  at  Stagnation  Point,  a  -  r$ 

0.053  in 

Assumed  Physical  Properties 

Sublimation  Temperature,  T 

Jfl 

1 200°R 

Diffusivity,  x 

1 .67  x  10  4  in^/sec 

Conductivity,  k 

3.61  x  10-4  B™ 

^  Ft  sec 

Density,  p 

74.4  lb/fr 

Specific  Heat,  c 

0.35  BTU/ 1 b°R 

Heat  of  Sublimation,  L 

450  BTU/ lb 
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Pre-Helt  Solution 

From  the  above  data  the  parameters  necessary  to  use  Fig.  2  can 
be  calculated.  These  are 


and 


0.569  • 


From  Fig.  2,  t^  =  3-30  x  10"5  , 


and  hence 

a2t' 

t  =  - IS.  -  2.96  seconds 

m  h 

A1 bat  ion  Solution 

The  ablation  parameters,  and  y,  calculated  from  the  given  data 
are 


<p  =  0.5H 

y  =  62.8 


The  non-dimensional  time  corresponding  to  the  conclusion  of  the  test  (t^) 
can  be  determined  from  the  given  data  and  from  the  pre-melt  solution.  That  is 


u? 
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n(t,-t  )  . 

-  •  m  _  0  -|Q  w  inmt* 

t  ^  j  3  •  78  1 0 


The  total  ablation  depth  predicted  by  the  theory  can  now  be 
determined  by  using  Fig.  3>  The  result  is 


a  -  rg  =  0.047  in. 

The  actual  ablation  depth  measured  at  the  stagnation  point  is  0.053  in. 
and  hence  the  theoretical  result  is  within  12%  of  the  experimental  result. 
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TIME  HISTORY  OF  ABLATING  SURFACE ,  ja-0.108 
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ABLATION  RATE  VERSUS  TIME ,  <fi  -  0.108 


